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Abstract. We derive the precise limit of SHS in the high activation en- 
ergy scaling suggested by B.J. Matkowksy-G.I. Sivashinsky in 1978 and by 
A. Bayliss-B.J. Matkowksy-A.P. Aldushin in 2002. In the time-increasing case 
the limit coincides with the Stefan problem for supercooled water with spatially 
inhomogeneous coefficients. In general it is a nonlinear forward-backward par- 
abolic equation with discontinuous hysteresis term. 

In the first part of our paper we give a complete characterization of the limit 
problem in the case of one space dimension. 

In the second part we construct in any finite dimension a rather large family 
of pulsating waves for the limit problem. 

In the third part, we prove that for constant coefficients the limit problem in 
any finite dimension does not admit non-trivial pulsating waves. 
The combination of all three parts strongly suggests a relation between the pul- 
sating waves constructed in the present paper and the numerically observed 
pulsating waves for finite activation energy in dimension n > 1 and therefore 
provides a possible and surprising explanation for the phenomena observed. 
All techniques in the present paper (with the exception of the remark in 
the Appendix) belong to the category far-from-equilibrium-analysis/far-from- 
bifurcation-point- analysis. 



1. Introduction 

The system 

, , dtu ~ Au = vf{u) 

dtv = -vf{u) , 

where u is the normaUzed temperature, v is the normahzed concentration of the 
reactant and the non-negative nonhnearity / describes the reaction kinetics, is a 
simple but widely used model for solid combustion (i.e. the case of the Lewis num- 
ber being +oo). In particular it is being used to model the industrial process of 
Self-propagating High temperature Synthesis (SHS). In the case of high activation 
energy interesting phenomena like instability of planar waves, fingering and screw- 
like spin combustion waves are observed. 
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In [2^ B.J. Matkowsky-G.I. Sivashinsky derived from a special case of ((ij a formal 
singular limit containing a jump condition for the temperature on the interface. 
Later it has been argued that the problem is for high activation energy related 
to a Stefan problem describing the freezing of supercooled water (see [20], O p. 
57]). Subsequently the Stefan problem for supercooled water became the basis 
for numerous papers focusing on stability analysis of ([T]), fingering, screw-like spin 
combustion waves etc. (see for example [9].|10j.[T2 ] .[llj.[T3].[8].[2]). 
Surprisingly there are few mathematical results on the subject: In [19| E. Logak-V. 
Loubeau proved existence of a planar wave in one space dimension and gave a rig- 
orous proof for convergence as the activation energy goes to infinity. 
Instability of the planar wave for a special linearization (and high activation en- 
ergy) is due to [4]. 

The present paper consists of three parts: in the first part we prove rigorously that 
in the case of one space dimension the SHS system converges to the irreversible 
Stefan problem for supercooled water (cf. Theorem 15. 1[) . In the time- increasing 
case we obtain also convergence in higher dimensions (see Theorem 14. ip . As the 
initial data of the reactant concentration enter the equation as the activation en- 
ergy goes to infinity, our result also seems to provide a possible explanation for the 
numerically observed pulsating waves (cf. [2], [H] and [5]). 

As a matter of fact, in the second part of our paper (Theorem 17. ip we use the 
spatially inhomogeneous coefficients in order to construct a pulsating wave for each 
periodic function u° (or Y^, respectively) on R", using the approach in [3]. We also 
obtain the spin combustion waves (or "helical waves" ) on the cylinder mantle (see 
Remark [7!2l). 

In contrast, we show in the third part (see Theorcm lS.ip that for constant v'^ in any 
finite dimension, no non-trivial pulsating waves exist. In addition, the formal result 
in the Appendix suggests that in one space dimension the planar wave is stable. 
Taken together (cf. section |9]) , our results strongly suggest a relation between the 
pulsating waves constructed in the present paper and the numerically observed 
pulsating waves for finite activation energy in dimension n > 1 (cf. [14j . [22] and 
[2]) and therefore provide a possible and surprising explanation for the phenomena 
observed. 

In the original setting by B.J. Matkowsky-G.I. Sivashinsky [201 equation (2)], ac- 
cording to our result Theorem 15. 11 



l,iV t oo the limit in this scaling is the solution of the heat equation; cf. Section 
Oand Theorem [53|) 




(3) 5tUoo - v°dtx = Awoo in (0, -f-oo) x f2. 



where are the initial data of Uoo and 
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In the SHS system with another scahng and a temperature threshold (see [2j p. 



= (1 - (TAr)iVyivexp((iV(l - (Tjv)(0jv - l))l{oN + (1 - (TN)eN))X{e^>e}, 
dtYN = -(1 - ctn)NYn cxp((Ar(l - cTN){eN - l))l{cJN + (1 - <yN)eN))X{e^>S} 

where 7V(1 — gn) » 1, (Jn G (0, 1) and 6 E (0, 1), each limit 6*00 of 6^ satisfies (cf. 
Section 16.21 and Theorem 15. ip 



To our knowledge this precise form of the limit problem, i.e. the equation with the 
discontinuous hysteresis term, has not been known. Even in the time-increasing 
case it does not coincide with the formal result in [20] ■ 

In the case that 9 00 (or Uoo, respectively) is increasing in time and (or K*^, re- 
spectively) is constant, our limit problem coincides with the Stefan problem for 
supercooled water, an extensively studied ill-posed problem (for a survey see [5]). 
As it is a forward-backward parabolic equation it is not clear whether one should 
expect uniqueness (see [6l Remark 7.2] for an example of non-uniqueness in a re- 
lated problem). 

On the positive side, much more is known about the Stefan problem for supercooled 
water than the SHS system, e.g. existence of a finger ([E]), instability of the finger 
([l8]). one-phase solutions ([6]) etc.; those results, when combined with our conver- 
gence result, suggest that similar properties should be true for the SHS system. 
It is interesting to observe that even in the time-increasing case our singular limit 
selects certain solutions of the Stefan problem for supercooled water. For example, 
u{t) ~ {k — l)X{t<i} + i^X{t>i} is for each k e (0, 1) a perfectly vahd solution of the 
Stefan problem for supercooled water, but, as easily verified, it cannot be obtained 



from the ODE 

1 /•* 

dtUs{t) = -dtexp{ — / cxp((l - l/(iis(s) + l))/£) ds) as e — » . 



Our approach does not involve stability or bifurcation analysis. For showing the 
convergence we use standard compactness and topological arguments. We prove 
that if the measure of the "burnt zone" is small enough in a parabolic cylinder, 
then in a cylinder of smaller radius there cannot be any burnt part. For the con- 
struction of pulsating waves we use the approach in [3] as well as blow-up arguments, 
Harnack inequality and so on. In order to obtain non-existence of pulsating waves 
for constant we use a Liouville technique. 

Let us conclude the introduction with a comparison to blow-up in scmilinear heat 
equations, as the main problem arising in our convergence proof, i.e. excluding 
"peaking of the solution" or burnt zones with very small measure, resembles the 
blow-up phenomena in semilinear heat equations. One could therefore hope to ap- 
ply methods used to exclude blow-up in low dimensions in order to exclude peaking. 



109-110]), 
(4) 



(5) dt9^ - y"9tx = in (0, -f 00) x ft, 

where Y'^ are the initial data of Fqo and 
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say in two dimensions. There are however problems: First, here, we are deahng 
not with a single solution but with the one-parameter family Ue concentrating at 
some "peak" as s gets smaller. Second, the e-problem is not a scalar equation but 
a degenerate system. Third, in contrast to blow-up, peaking would not necessarily 
imply going to +00. Fourth, our limit problem is a two-phase problem while 
most known results for blow-up in semilinear heat equations assume the solution 
to be non-negative. Fifth, in our problem it does not make much sense studying 
the onset of burning, say the first time when > — e, whereas studying the time 
of first blow-up can be very reasonable for semilinear heat equations. The last and 
most important difference is that while semilinear heat equations are parabolic and 
therefore well-posed in a sense, our limit problem contains a backward component 
making it ill-posed. 

2. Notation 

Throughout this article R" will be equipped with the Euclidean inner product 
X ■ y and the induced norm |a;| . Br{x) will denote the open n-dimensional ball of 
center x , radius r and volume r"u;„ . When the center is not specified, it is assumed 
to be 0. 

When considering a set ^ , xa shall stand for the characteristic function of A , 
while V shall typically denote the outward normal to a given boundary. We will 
use the distance pardist with respect to the parabolic metric d{{t, x), {s,y)) = 
v/|t-s| + |x-y|2. 

The operator dt will mean the partial derivative of a function in the time direction, 
A the Laplacian in the space variables and £" the n-dimensional Lebesgue measure. 
Finally W^' ^ denotes the parabolic Sobolev space as defined in [T7]. 

3. Preliminaries 
In what follows, is a bounded C^-domain in R" and 

Te(0,+oo) 

is a strong solution of the equation 

(Q\ dtUe{t, x) - Aue{t,x) = -v°{x)dt exp{-^ J* ge{us{s,x)) ds) , 

^ ' Ue(0, •) = u° in n, Vue • 1/ = on (0, +00) x dfl ; 

here (7e is a non-negative function on R satisfying: 

0) 5e is for each e € (0, 1) piecewise continuous with only one possible jump at zq, 
(7e(zo~) = Qeizo) = in case of a jump, and satisfies for each e € (0, 1) and for 
every z G R the bound ge{z) < Ce(l -I- |z|). 

1) ge/e as £ ^ on each compact subset of (—00, 0). 

2) for each compact subset K of (0, +00) there is ck > such that Tahi{g^, ck) — > ck 
uniformly on /\ as e — > 0. 

The initial data satisfy < < C < -foo, w° converges in L^{Vl) to as £ ^ 0, 
(u^)Eg(o,i) is bounded in it is uniformly bounded from below by a constant 

Mmin, and it converges in L^{^l) to u° as e 0. 

Remark 3.1. Assumption 0) guarantees existence of a global strong solution for 
each e e (0, 1). 
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where x{t, x) 



4. The High Activation Energy Limit 

The following theorem has been proved in [21j . Let us repeat the statements 
and its proof for the sake of completeness. 

Theorem 4.1. The family (iie)e£(o.i) is for each T G (0, +oo) precompact in 
L^((0,T) X fl), and each limit u of (uj)£g(o.i) '^■s ^ sequence Em ^ 0, satisfies 
in the sense of distributions the initial-boundary value problem 

(7) dtu - v^dtx = Au in (0, +oo) x ft, 

u{0, ■) = u° + v°H{u°) mn , Wu ■ ^ on {0, +oo) x dn , 

e[0,l], esssup (o^t)u(-,a::) < , 
= 1, esssup (Q > , 

and H is the maximal monotone graph 

r = 0, z < 0, 
Hiz)\ G[0,1], z = 0, 
[ = 1, z > . 

Moreover, x is increasing in time and u is a supercaloric function. 

If (we)eg(o,i) satisfies dtu^ > m (0,r) x 17, then u is a solution of the Stefan 

problem for supercooled water, i.e. 

dtu — v^dtH{u) = Au in (0, +oo) x . 

Remark 4.2. Note that assumption 1) is only needed to prove the second statement 
"//....". 

Proof. Step (Uniform Bound from below): 

Since is supercaloric, it is bounded from below by the constant Umin- 
Step 1 (L2((o,T) X r2)-Bound): 

The time-integrated function Vi,{t,x) := j^Us{s,x) ds, satisfies 

(8) dtvS, - Av,{t, x) = w,{t, x) + u°{x) 

where is a measurable function satisfying < < C*. Consequently 

Jn ^ Jn Jo Jn 

^ Jo Jn ^ Jn 

implying 

i-T I- r 

,0h2 



(9) 



/ j ul< Thc+\ul\ 
Jo Jn Jn 



Step 2 (L2((o,T) X f})-Bound for Vmui{u^,M): 
For 

and any A/ e N, 



Mz - AP/2, z>M 



Gm{u,)^Gm{u",) + [ / |Vmin(ii„Af)|2 
Jo Jn 
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1 -v"mm{us,M)dtexp{~- / g^{ue{s,x)) ds) . 
Jn ^ JO 

As dt exp(— i g^{ue{s, x)) ds) < 0, wc know that dt exp(— i J* g^{u^{s, x)) ds) is 

bounded in i°°(r2;Li((0,T))), and 

T r I rt 

-w^ niin(we, Af)9t exp( / g^{u^{s, x)) ds) 



Jn 



e 



< C sup max(min(ue,M),0) < CMC^n) . 
Jn (o,T) 

Step 3 (Compactness): 

Let xm : R. ^ R be a smooth non- increasing function satisfying X(_oo,a/-i) ^ 
Xm < X(-oo,A/) ^J^id let (f>jv/ be the primitive such that $m(-z) = z for z < M ~ 
1 and $A/ < -A'^- Moreover, let {(l)s)se{o,i) be a family of moUifiers, i.e. (ps G 
Co'^(R"; [0, +oo)) such that J 4>5 = I and supp 0^ C ^^(O) . Then, if we extend 
and by the value to the whole of (0, +oo) x R", wc obtain by the homogeneous 
Neumann data of that 

dti<^M{Ue)*(l)s){t,x) 



XM{ue) ^XnAwe - v°dtexp{~^ J 



XM{u^){t, y) xn{y)^u^{t, y) 



1 



v°{y)dtexp{-- J ge{u^{s,y)) ds)j(j>s{x - y) dy 
(l>s{x x'Mi'Ue{t,y))xn{y)\"^u,{t, y)p 



1 



Consequently 



/R 

- XAf(?ie(i,y))w"(y)9iexp(-- / g^{u^{s,y)) ds) 

£ Jo 

+ XM{ue{t,y))xn{y)^u^{t,y) ■ W(j)s{x - y) dy . 



T 

|5t($MK)*05)| < Ci(17,C,M,5,T) 



JR" 

T 



and ^ 

/ / |V($A/(Me)*05)| < C2(r!,A/,^,T) . 

Jo JR" 

It follows that {^M{ue) * 4>s)ee{o 1) is for each {M,S,T) prccompact in L^{{Q,T) x 
R"). 

On the other hand 

1 

T r ( i-T 

2 



\^m{u,) * 05 - $a/(h,)| < C3 / |V$A/(^^e)r 

^R" \ Jo Jn J 

+ 2{M - u„,in)r/:"(B5(9r!)) < Ci{C,^l,u^i^,M,T) 6 . 
Combining this estimate with the precompactness of ($Af (ug) * 4>s)e£(o.i) obtain 
that <i>Af (we) is for each (M, T) precompact in L^((0, T) x R"). Thus, by a diagonal 
sequence argument, we may take a sequence ^ such that $Af(we„) zm a.e. 
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in (0, +00) X R" as m 00, for every M £ N. At a.e. point of the set {zm < M — 
1}, Ue^^ converges to zm- At each point {t, x) of the remainder PlA/eNi-^A^ — ^'^^ 
the value u^^ {t, x) must for large m (depending on (M, t, x)) be larger than M — 2. 
But that means that on the set CluiGNi^^^ > M — 1}, the sequence (u£^)„jg]sj 
converges a.e. to +00. It follows that (ue^)„jg]\f converges a.e. in (0, +00) x 
to a function z : (0, +cx3) x 17 — > R U {+cx)}. But then, as is for each 

T S (0, +00) bounded in L^((0,T) x f2), (uj„J,^g]vf converges by Vitali's theorem 
(stating that a.e. convergence and a non-concentration condition in imply in 
bounded domains L^'-convergence) for each p G [1, 2) in L^((0, T) x fJ) to the weak 
L2-limit u of (u£,J„eN- I* follows that C^+HClMe-Ni^M > M-1}) = C"+^{{u = 
+00}) = 0. 

Step 4 (Identification of the Limit Equation in esssup (p t)^ > 0): 

Let us consider (t, x) E (0, +00) x il such that u^^ (s, x) — > u{s, x) for a.e. s e (0, t) 

and u{-,x) £ i^((0, i)). In the case esssup (Q_i)u(-,a;) > 0, we obtain by Egorov's 

theorem and assumption 2) that exp(— g^,^ {ue„^ {s, x)) ds) ^ as m — > 00. 
Step 5 (The case dtu. > 0): 

Let {t,x) be such that u^^(t,x) —^ u(t,x) = A < 0: Then by assumption 1), 



exp( / g^^ [Ue^ [s, X)) as) > exp(— t ■ — ) ^ 1 as to ^ cx) . 

Jo ^rn 

Remark 4.3. We also obtain a rigorous convergence result in the case of (higher 
dimensional) traveling waves with suitable conditions at infinity. In this case our 
(W^'^)-estimate (Step 2) implies a no- concentration property of the time- derivative. 

5. Complete characterization of the limit equation in the case of 

ONE space dimension 

The aim of this main section is the following theorem: 

Theorem 5.1. Suppose in addition to the assumptions at the beginning of Section 
that the space dimension n = 1 and that the initial data converge in to 
a function satisfying Vu" ^ on {vP = 0}. Then the family (Me)eg(o.i) is for 
each T G (0, +00) precompact in L^{{0,T) x Q,), and each limit u of (ue)£g(o.i) 
a sequence Em ^ 0, satisfies in the sense of distributions the initial-boundary value 
problem 

(10) dtu - v^dtx = Au in (0, +00) x 17, 

it(0, •) = M° + v"H{u") m , Vu • 1/ = on (0, +00) x dil , 
where H is the maximal monotone graph 

f = 0, z < 0, 
H{z) I e [0,1], z = 0, 
[ = 1, z > 

{=0, esssup (o j)u(-,a;) < , 
e[0, 1], esssup (o.t)w(': 2;) = , 
= 1, esssup jju(-,a;) > , 

Although we assume the space dimension from now on to be 1, we keep the 
multi-dimensional notation for the sake of convenience. Moreover we extend by 
even reflection at the lateral boundary to a space-periodic solution on [0, -l-oo) x R. 
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X 



Figure 1 . Clearing out 
We start out with some elementary lemmata: 

Lemma 5.2 (Clearing out). There exists a continuous increasing function to : 
[0, 1) [0, +oo) such that uj(0) = and the following holds: suppose that k < 0, 
that e < uj{k), that 6 £ (0, 1) and that < (1 + lu{S))k on the parabolic boundary 
of the domain Q{to, S, 01, 4>2) ■= {(t, x) : < to — 26 < t < Iq, (f>i(t) < x < 4>2{t)}, 
where (j>i < 02 are C'^ -functions. Then < k in Q{tQ,S,(j>i,(j)2) (cf. Figure\^. 

Proof. Comparing in Q(to, (5, 0i, 02) to the solution of the ODE 

V'{t) = CgM/e. y{to - 26) = (1 + LoiS))K 
we obtain the statement of the lemma. 

Lemma 5.3. For almost all k < the level set ((0,+oo) x il) n {u^ — k} is a 
locally finite union of -curves. For such k we define the set 

S'„,£ := {{t,x) e (0,+cx)) X n) : Ue{t,x) > k 
and there is no (to, 6, 0i, 02) S [0, +oo) x (0, 1) x x 
such that < K on the parabolic boundary of the domain Q{to, i5, 0i, 02)} 

(cf. Figure\^. Then dS^.e = [j'j'jli graph {gj^K,^e) where gj^K,,e '■ [0,Tj,K,e] R are 
piecewise -functions and iV^.e is for small £ bounded by a constant depending 
only on the limit u° of the initial data. 

Remark 5.4. For illustration of the definition of S^^e, imagine the set {u^ > k} 
filled with water in a {t, x) -plane where t represents the height. Our modification of 
{u^ > k} means then that the water is now allowed to flow out through the "bottom" 
{t = 0}. 

Proof of Lemma \5.3l By the definition of S'^^e and by the fact that is su- 
percaloric, each connected component of dS^^e is a piecewise C^-curve and touches 
{t = 0}. Therefore the number of connected components is for small e > bounded 
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t 




X 



Figure 2. The set S, 



by a constant N depending only on the Umit vP of the initial data. 
Let us consider one connected component 7 of DSk^s- By the definition of S^.e and 
by the fact that is supercaloric, the derivative of the time-component of the 
piecewise C^-curve 7 can change its sign at most once! Thus we can define for each 
curve 7 one or two piecewise -functions of time such that 7 is the union of the 
graphs of the two functions. The total number of graphs A^^.e is therefore bounded 
by 27V. 



Proof of Theorem \5.U 

By Theorem l4.1l we only have to prove that x = in the set {esssup (o,t)w(-, x) < 0}. 
The main problem is to exclude "peaking" of the solution Ug, i.e. tiny sets where 
> n. Here we show that in the case of one space dimension, "peaking" is not 
possible. More precisely, if the measure of the set > k is small in a parabolic 
cube, then Ue is strictly negative in the cube of half the radius, uniformly in e. The 
proof is carried out in two steps: 

Step 1: Let (£m)„g]Sf be the subsequence in the proof of Theorem 14.11 As a.e. 
point (t, x) S ((0, -l-oo) X R) n {u < 0} is a Lebesgue point of the set {u < 0}, we 
may assume that there exists k < such that for any 6 e (0, 1), sufficiently small 
ro > and every e„i G (0, Eq), 

£^{{{t~2ro,t) X B2roix)) n {u,^ < 2k}) > 9C^{{{t - 2ro,t) x B2,,,(a;))) . 

Step 2: Suppose now that ((t — ro, i) x Br.g{x))r\{ue^ > «;} ^ (where k is chosen 
such that = k} and = 2k} are locally finite unions of C-'^-curves): then 

((t — ro, t) X Brg{x))r]dSK,,e„, and ((t — 7'o, t) x Bro{x))r\dS2K,e„^ must by Lemma [Ol 
be connected to the parabolic boundary of {t — 2ro , t) x B2ro (x) ■ The ( (0, T) x fi)- 
Bound for Vmin(Me,M), the fact that ^^(((t „ 2ro, t) x ^2^0 (a:)) n {u^^ < 2k}) > 
6'£^(((i— 2ro, t)xB2ro{x))) and Lemma [5?2l implv now (see Figure[3]) that there must 
be an "almost horizontal" component of 88^,5^ (cf- Figure [3]) with the following 
properties: 

for any S £ (0, 1), there are t — tq < ti < t2 < < t such that (see Figure [5]) 
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t , 














r 

J(t-ro,t) X B,.g(x) 











X 



Figure 3. Situation excluded by the L^(VF^'^)-estimate 




> K 



Figure 4. The main task is to exclude almost horizontal propagation 

— ti ^ a.s Em ^ 0, for some j 

|5i,K,E„(i2) - > ci > , 

and 

^\{y e Br.ix) ■■ ue,jh,y) > 24) < s , 

/ \u,Jh,y)\dy <S. 

JBr„ (a;)n{«e„,(t3,-)>2«;} 

We may assume that ci < tq, that gj.K.e^ih) = sup^^^^^^^ 5j,k,£™: that g.j^K.,e^{h) > 
5i,K,em(^i) siiid that Ue„^{s,y) > k for some d> and {s,y) G (^1,^2) x Brg{x) such 
that gj,K.£„(s) <y < + .gj^K.e,„(s). We define the set D^^^ := 

{(s,?/) : ii < s < ta,?/ < 9j,K,e„Ss) for s G (^1,^2) and y < gj^^.c^ih) for s € [<2,i3)} 

(cf. Figure[5]) and the cut-off function := max(0,min(?/— gj^K,e^(t2) — 
y)). It follows that 

cIk/A + 2S + o{1) > 0(1) + V(y)(we„(i3,2/)-«:)rfy 
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Figure 5. The set D^^ 



/■Sj.K.em (*2) rt2 

> o{l)+ (j){y)ue^{t3,y) dy - K '/'(5j,K,e„. (s))5j,^_^„, (s) ds 

•^Sj.i^.em (*l) •'tl 
/■Si.K.Em (*2) /"'S 
Sj.K.em (*l) 



JDc.^ -Id, 



pt2 



> - \7ue„^ as ^ , 

a contradiction for small provided that S has been chosen small enough; in the 
third inequality we used Lemma [5T^ and the convergence to is due to the uniform 
L2(VFi.2)_bound. 



6. Applications 

Although the limit equation is an ill-posed problem, the convergence to the limit 
seems to be robust with respect to perturbations of the e-system and the scaling: 
here we mention two examples of different systems leading to the same limit. Other 
examples can be found in mathematical biology (see [H] and [M])- 
For the convergence results below we assume that the space dimension is 1. 
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6.1. The Matkowsky-Sivashinsky scaling. We apply our result to the scaling 
in [201 equation (2)], i.e. 

. dtujy ~ Aujv = (1 - crN)NvN cxp(7V(l - l/u^v)), 

^ ' dtVN ^~Nvmcxp{N{1-1/un)), 

where the normalized temperature un and the normalized concentration vn are 
non-negative, (o'Ar)^g]\g- CC [0, 1) (in the case (Jat t l,iV "f oo the limit equation in 
the scaling as it is would be the heat equation, but we could still apply our result 
to tiAr/(l — (Tn)) and the activation energy N ^ oo. 
Setting Umin := ^/N,Ue := un — 1 and 

r exp((l-l/(z + !))/£), z>-l 
' \ 0,z < -1 

and integrating the equation for vn in time, we see that the assumptions of Theorem 
15.11 are satisfied and we obtain that each limit Uoc, cToo of uat, ctat satisfies 

(12) dtUoQ - (1 - croo)^'°9*-H^(esssup (o^()Uoo) = Amqo in (0, +oo) x VL, 

Uoo(0, •) = + v°H{u°) in n , Vuoo • = on (0, +oo) x dn , 

where iP are the initial data of Voo- Moreover, x is increasing in time and Uqq is a 
supercaloric function. 

6.2. SHS in another scaling with temperature threshold. Here we consider 
(cf. m p. 109-110]), i.e. 

(13) 

dtdN — ^(^N 

= (1 - aN)NYN exp((7V(l - aN^ON ~ l))/{aN + (1 - <TN)0N))X{e«>S}, 
dtYN - -(1 - aN)NYN exp((iV(l - ctn^On - l))/i<JN + (1 - aAr)0jv))X{e„>e} 

where A^(l — ctat) >> l,a]s[ G (0,1) and the constant 6 E (0,1) is a threshold 
parameter at which the reaction sets in. 

Setting Uinin = -1,^ := l/(iV(l - (Tn)), ^(e) := 1 - aN, u^ := On - 1, 

„ / exp((z/(K(£)z + l))/e), z > ^ - 1 

^"^^ \ 0,z<^-l 

and integrating the equation for in time, we see that the assumptions of Theorem 
15.11 are satisfied and we obtain that each limit Uoc of un satisfies 

(14) dtUao - v°dtH{csssiip (o,t)Uoo) = Auoo in (0, +oo) x il, 

Uoo(0, ■)^u" + v°H{u") in , Vuoo • ^ = on (0, +oo) x dn , 

where v'^ are the initial data of Woo- Moreover, x is increasing in time and 
supercaloric function. 

7. Existence of Pulsating Waves 

The aim of this section is to construct pulsating waves for the limit problem. 
For the sake of clarity we have chosen not to present the most general result in the 
following theorem. Moreover we confine ourselves to the one-phase case. 
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Theorem 7.1. (Existence of pulsating waves) 

Let us consider a Holder continuous function defined on R" that satisfies 

> 1 and v° {x + k) ^ v° (x) for every k e Z'\ .x G R" . 

Given a unit vector e € R" and a velocity c > 0, there exists a solution u{t, x) of 
the one-phase problem 



(15) 

which satisfies 
(16) 



dfU > and — ^lq — Jjp < u < , 



u{t,x + k) = u{t - ^—^,x) for every fc e Z", (t, x) G R x R" 



u{t, x) = for X ■ e — ct < and limsup u{t, x) = —fio i 

x-e~ct — >+oo 

where the last limit is uniform as x ■ e — ct tends to +oo . 

Remark 7.2. By modifications of the following proofs and of the theory in [3], 
it is possible to replace R" in Theorem \7.1\ by a smooth source manifold. Taking 
for example x R, we obtain the screw-like pulsating waves observed in spin 
combustion (also called "helical waves"; see for example |14j . [22j . [2].[14j J. 

Let us transform the problem by the so-called Duvaut transform (see [13]), set- 
ting w{t,x) = — J^°°u{s,x) ds. In this section we will prove the existence of a 
pulsating wave w. More precisely, Theorem 17.11 is a corollary of of the following 
result which will be proved later. 

Theorem 7.3. (Pulsating waves for the obstacle problem) 

Under the assumptions of Theorem \7.1\ there exists a function w{t,x) solving the 
obstacle problem 

dtw = Aw — X{iu>o} on R X R" , 



(17) 

with the conditions 
(18) 



w > , —fj,o < dtw < 0, dttw > 



w(t,x + k) = w(t - ^—^,x) for every fc G Z", (t,a;)GRxR", 
c 

w{t, x) = for X ■ e — ct < and dtw{t, x) — > — /io as x ■ e — ct — > -|-oo . 
The convergence is uniform as x ■ e — ct tends to -\-oo. 
Proof of Theorem I7.il 

Simply set u(i, a;) := dtw{t, x) with w given by Theorem l7.31 and use the fact that 
X{ii<o} = X{w>o}- To check this last property, it is sufficient to exclude the case 
where w{to,xo) > and dtw{to,xo) = at some point {to,Xo)' using the fact that 
dfW is caloric in {w > 0} as well as the strong maximum principle, we deduce that 
dtw{t,x) — for t G (— oo,<o] said a; in a neighborhood of xq. This contradicts the 
last hne of (UHl). 
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Proof of Theorem [73 

We will prove the existence of an unbounded solution w in six steps, approximat- 
ing w by bounded solutions of a truncated equation, for which we can apply the 
existence of pulsating fronts due to Berestycki and Hamcl [3] . 

Step 1: Approximation by bounded solutions and estimates of the ve- 
locity 

For any < ^ < M, let us start by approximating the function X(o,+oo) by the char- 
acteristic function g = X{o.A]- In that case we can compute explicitly the traveling 
wave (0, Co) (unique up to translations of 0) of 

(19) CQ(f>' = (f)" - g{(f)), 0' < on R, ^(-oo) = M and (j){+oo) = : 
Let us define for co > 0, M > 0, sq G (— oo, 0) and si G (so, 0) 

M (1 - e=o(^-'^")) for s E (-oo, si] , 

0(s) ^ <( \ (e'=°" - 1 - Cos) for s G [si, 0] , 

] for s G [0, +oo) . 

For any A G (0, M) and for suitable so, si, we see that (j) is continuous and satisfies 
0(si) = A, which fixes the parameter si as a function of A. Moreover we see that 
4> is of class if and only if si = — coM and 

(20) M-A^\(l- e-^"^^ 

Cq \ 

Thus A is determined in terms of the velocity co and M. The above calculations 
show in particular that (l>{cot—e-x) is a good bounded approximation of the solution 
of P?)) - pB]) in the case = 1, i.e. the travehng wave case. 
In all that follows let A be given by (PO)) . 

Now, when the function g in (|19[) is replaced by a Lipschitz continuous function 
whose support is a compact interval, there are known results on the existence of 
pulsating waves. For such g, it is possible to apply Theorem 1.13 of Berestycki, 
Hamel , which states the existence (and uniqueness up to translation in time) of 
bounded pulsating solutions traveling at a unique velocity. Bearing that in mind, 
we define gM as a Lipschitz rcgularization of the characteristic function g such that 
supp gM = [0, ^] and - for later use - 

9M = l on [l/M,A/2], 0<gM<l on R, 
^ ' and .g;, >0, .g^, < on (0,^/2). 

Let us call c^j the unique velocity of the traveling wave equation (|19p with g replaced 
by gM- As (0, co) can be shown to be unique up to translations of 0), c^^ — ^ co as 
gM 9- 

Then there exists by Theorem 1.13 of |3] a bounded pulsating wave wm traveling 
at velocity cm such that 

dtWM = AifA/ - v^gM^WM) on R x R" , 

dtWM < and limsup WM{t,x) = < wm < M = liminf WM(t,x) 

x-e—CMt — ► — C50 x-e—CMi *-t-oo 

and WM{t,x + k) ^ WMit ~ ^—^,x) for every fc G Z", (t,x)GRxR". 

Cm 
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From the assumption 1 < v'^ and the comparison principle (see Lemma 3.2 and 3.4 
of [3]) we infer that the velocities satisfy the ordering Cj[j < cm. Similarly, defining 
A = IIw^IIl^ and comparing to WM{^t, V^x), we get cm < c°/"\/A. 
Furthermore the above comparison principles tell us that the velocity cm (resp. 
c^/) is continuous and non-decreasing in AI. 

For all that follows let c > be an arbitrary but fixed velocity for which wc want 
to construct the pulsating wave. Then, for any M > we can adjust A e (0, AI) 
such that for cq defined by (PO)) . 



c = CM e [co/2, 2coVII«"IIl-] • 

In order to pass to the limit as M +oo, we need to get some bounds on the 
solution first. To this end, rotating space-time proves to be very convenient: 
Step 2: Space-time transformation and first estimates on the time deriva- 
tives 

Let us introduce the function ium defined by WMis^x) = wm C^'^^'^ , x) which is 
periodic in x and satisfies 

LwM ^ v°{x)gM{wM) with Lwm ^ ^wm + dssWM " 2de.sWM " cd^WM 

and linis^-oo wm(s, x) — AI, linis^+oo u'a/(s, a;) = uniformly with respect to 
X. Using pTj) . we obtain LdgWM — cidgWM = and LdssWM — cidssWM < on 
{wm < A/2} for ci{s,x) = v°{x)g'^j{wM{s,x)) > on this set. We deduce from 
the maximum principle (see Lemma 3.2 and 3.4 of [3]) that for any sq G R such 

that supj^^ _^^)^p^. WM < A/2, 

inf[,o,+oo)xR" 9sWM = infR" dsWM{so, •) 

and infj^^^^^-i^jj^r. 9ssWm = min (0, infjj^-. 5ssWa/(so, •)) ■ 

Step 3: Bound of the solution from above 

From the fact that gM is bounded by 1 , and from the Harnack inequality, we deduce 
that there exists a constant Ch G (1i+oo) such that for any r > and for any 
point {to,xo) 

(23) sup WM{to-r'^,-) <Ch { hif WM{to, ■) + r"^ X 

where A = For wm that means that - setting so = do — e ■ xo - 

sup WM{so-cr'^-e-y,xo+y) < Ch [ inf wm(so - e • xo + y) + r^A 
yes^/2(o) \yeB^/2(o) J 

for r > ^/n/2. We will now use the fact that the unit cell (—1/2, 1/2)" is contained 
in the ball 5^/2(0). Using first the monotonicity of wi\j in the variable s, and 
second the periodicity of wm{t, y) in y, we get for tq sq — \fn/2 

(24) supiyM(To - cr^ ^fh, ■)<Ch\ inf wm{tq, •) + r^A ) . 
R" VR" / 

By a translation in time wc may assume that 

(25) = inf{T: wm{.s,x) <\/AI for s > r, a; G R"} 
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and get the bound 

(26) wj\/(s, x) < max (1/Af, a — /3s) 

for some constants a,/3 E (0, +oo) and every large positive M. 

Step 4: Passing to the limit 

By estimate ([25)) . we can pass to the hmit as M +oo and obtain AI — A ^ I/cq- 
Moreover, passing to a subsequence if necessary, wm converges locally in W^'^ to 
li satisfying 

dtw < and limsup w{t,x) ~ < w 

x-e—ct — ' — oo 

and w(t, X + k) ^ wit - - — , x) for every k € Z", (i, a;) G R x R" . 

c 

Furthermore, we obtain for w related to w by w{s,x) = w ( "^^^'^ , x) that 

wt = Aw - v^X{w>Q} ; 

here we used the fact that w, being locally a Wp'^-function, satisfies dtW = = Aw 
a.e. on the set {w = 0}. 

In order to conclude dttw > the following non-degeneracy property will prove 
to be necessary: 

Step 5: Non-degeneracy property and bound from below 

Let us assume that WM{to, xo) E {l/AI, A/2). Using the fact that v'^{x)gM{z) > 1 
for z g [l/M,A/2\, we can use the usual parabolic maximum principle, comparing 
max(wM, 1/-^) to the function 

h{t,x) = WM{to,xo) + -^[x - + -j-ito - t) 

on the set 

{1/M < WM < A/2} n Q- {ta, xq) , 

where Q~{to,xo) = {{t,x) : to — < t < to, \x — xo\ < ?'}. We get for every 
r > the following non-degeneracy property: 

(27) sup WM > min (wM{to,xo) + -^r^, A/2 

Combined with the Harnack-type inequality (|24p for some radius r' < r, we obtain 
the following bound from below: 

(28) wm(s, x)>a'>0 for s < si < 
for some constants a' and si and every large M . 

Step 6: Further estimates on the time derivative of the limit solution 

By the bound from above in Step 3, we obtain that 

(29) \w{t,x)\<Ci + C2m + \x\) , 

where Ci and C2 are finite positive constants. Let now {tk,Xk) S {w > 0} be a 
sequence such that 

ctk - Xk ■ e ^ -00 . 

Then by the result in Step 5, 



(30) dk := pardist ((ifc, Xfc), d{w > 0}) > C3 vW+H^ 
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for some constant C3 > 0. So if is a solution of dtw — Aw = — in QdkitkjXk)- 
Defining 

w{tk + dlt, Xk + dkx) 
Zk{t,x) := -2 , 

(f29l) and (|30|) imply that Zk is a solution of dtZk — Azk = —v^{xk + dkx) in Qi(0) 
satisfying 

sup \zk\ < C4 , 
Qi(o) 

where C4 is a constant not depending on k. Consequently dt'w{tk,Xk) — dtZk(0) is 
bounded, implying that limsup(t_^)g{^>o} \dtw{t,x)\ < +00. Passing if 

necessary to a subsequence, we obtain by the periodicity of v'^ a limit z satisfying 
dtz — Az ~ — /io in Qi(0). Moreover we infer from the fact that w is periodic in the 
space variables that z is constant in the space variables. Thus dtz = — /ig in Qi(0). 
From regularity theory of caloric functions it follows that 

lim dttw ~ . 

But then a combination of the comparison principle (|22p and of ((^51) yield 

-Mo <dtw <0 on R X R" 

and 

dttW > on R X R" . 
This ends the proof of the Theorem 17.31 

8. Non-existence of pulsating waves in the case of constant initial 

concentration 

We consider solutions u of the one-phase limit problem with constant initial 
concentration in any finite dimension, i.e. 

(31) dtu ~ dtX{u>o} = Am in R x R" , 

and prove that u cannot be a non-trivial pulsating wave in the sense of (lisp . ()16p . 
More precisely: 

Theorem 8.1. (Non-existence of pulsating waves for constant initial con- 
centration) 

Let u be a solution of US]). U6\) in dimension n > 1 with v*^ = constant. Then 
u(t, x) = u{t — e ■ x/c, 0), i.e. u is a planar wave. 

Proof of Theorem \8.1\ 

We set w{t, x) = — J^°° u{s, x) ds > 0. From the proof of Theorem 17.11 we know 
that w > a and only if u < 0. As dtu > we obtain 

dtw = Aw- v°X{w>o} , 
and w satisfies ([T7)) . For any ^ G R", we define the "tangential difference" 

z^{t, x) = w{t , X — £,) — w{t, x) 

c 

which satisfies 

(32) {dt - A)z^ ^ -az^, 
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if z^it,x) = 

V Ui(i-5i.l-{)-tl,(i,l) y 
From PT)) and the definition of wc infer that 

Iftz^l < 2^*0 = 2w° in R"+i , 

dtz^{t, x) — > uniformly in t, a:, ^ as ct — e ■ x — > — oo . 

Moreover (fT8|) and (|28p as well as the definition of tell us that for some sq g 
(0, +oo) not depending on ^, 

(at-A)z« = in e • x| > So} . 

Furthermore we obtain from the comparison principle (see Lemmata 3.2 and 3.4 in 
[3|) that 

(33) |atz«(t,a;)| < 2w°e'=*-"-=^+"'' 

and — integrating this estimate for t e (— oo, £a±£:£) and using that = in 
t > ''o+J''^ — we obtain that z^ is bounded on R"^^ by a constant not depending 
on ^. 

Liouville's theorem for the heat equation implies therefore that for each sequence 
{tm,Xjn) such that ct + •) converges locally uniformly 

in R""^^ (and uniformly with respect to ^) to a constant K depending on the choice 
of ^ and the sequence {tm,Xm)- As we know that /^_|_[g z^{t, y) dy = for every 
(t, x) e R"+^ (sec (dH)), it follows that K = and that 

z^{t + ■,x + ■) ^ locally uniformly in R"'^"'^ as ct ~ e ■ x — oo ; 

the convergence is also uniform with respect to ^. 
Finally we define 

r]{t,x) := sup |2:^(t,a;)| . 
?gR" 

The function t] is by (|32|) a bounded subcaloric function. Moreover, by construction, 

dyT]{t - ^-—,2; - y) = . 

But then rj{t,x) = f{ct - e • x), cf - /" < in R, / € Wl;^CR), lim,^_oo f{s) = 
lim^^+oo /(s) = and / is bounded from above, implying that / = 0, that = 
and that w{t — x — ^) — w{t, x) = for every t € R and x, ^ G R". We obtain 
the corresponding result for u. 



9. Conclusions 

In this section we try to take a conclusion of the previous sections. 
Let us consider a sequence of solutions (ue,!)^) of the e- problem ^ satisfying for 
example the assumptions in the time-increasing case of Theorem I4.H and suppose 
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that they are getting closer and closer to one-phase pulsating waves as e — > 0, i.e. 

for some —>■ +00, 

(34) 

{e ■ k 
uJt, x^k) = oil) + Ue{t , x) for every k G Z", {t,x) G {t^, +00) x R' 
Ui,{t, x) = 0(1) for X ■ e ~ ct < —t^ and lim limsup ite(i, x) = — /io • 

Translating to 0, wc obtain by Thcorem l4.1l a sequence (^^(tg -\-t,x), v^{ti; -\-t, x)) 
that is locally compact in and converges to u satisfying all assumptions in 
Theorem 17. II except the assumptions for v^{x). Let us assume t;*^ > c > 0, i.e. there 
has been enough fuel everywhere initially, and let us show that is Z"-periodic: 
for the integrated function w as in section [7] and k G Z", 

= (dt- A)(w(t,x + k) -w{t-^-^,x)) 

c 

= v°{x)X{w{t-^,x)>o}~'"'^i.^ + k)x{w(t,x+k)>o} = {v°{x)-v°{x + k))x{w(t,x+k)>o} 
and therefore, choosing —t large, v^{x) — v^{x + /c) = 0. 

But then Theorem 18.11 tells us that u must either be one of the pulsating waves 
constructed in section [7l or a planar wave. 

So our results when combined, strongly suggest a relation of the numerically ob- 
served pulsating waves for dimensions n > \ and the pulsating waves constructed 
in section [71 Of course there remain possible alternatives, for examples our re- 
sults do not exclude the non-existence of pulsating waves for small e, or the strong 
convergence of pulsating waves to planar waves from far away, or the existence of 
non-trivial true two-phase pulsating waves in the case of constant We leave 
these possibilities open to mathematical discussion. However in communication 
with mathematicians doing numerics for pulsating waves of the very systems men- 
tioned in the introduction the above alternatives have been considered unlikely. 
Let us also mention that the chaos reported in [l] for finite activation energy 
presents no contradiction to our results: we never claimed that the dynamics around 
the pulsating waves is simple. 

On the other hand, our result - proving rigorously existence of pulsating waves 
for spatially inhomogeneous coefficients and non-existence for constant coefficients 
in any dimension - suggests that numerics - introducing nolens volens spatially 
inhomogeneous coefficients - makes the limit problem more unstable (see also the 
Appendix), and that grid information etc. will strongly influence the numerically 
obtained pulsating waves unless appropriate counter-measures are taken. 

10. Open questions 

The most pressing task is of course to study for space dimension n > 2 the 
existence or non-existence of "peaking" of the solution in the negative phase. A 
related question is whether {ug)^fz^o,i) is bounded in L°° in the case of uniformly 
bounded initial data. Although this seems obvious, it is not clear how to prevent 
concentration close to the interface. 

Another challenge is to use the information on the limit problem gained in the 
present paper to construct pulsating waves for the e-problem. 
Uniqueness for the limit problem (the irreversible Stefan problem for supercooled 
water) in general seems unlikely. One might however ask whether time-global 
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uniqueness holds in the case that u is strietly increasing in the xi-direction. By the 
result in [7] for the ill-posed Hele-Shaw problem, time-local uniqueness is likely to 
be true here, too. 



11. Appendix: Formal stability in the case of one space dimension 

AND constant INITIAL CONCENTRATION 

We consider solutions u of the one-phase limit problem with constant initial 
concentration in one space dimension, i.e. 

(35) dtu - dtX{u>o} = dxxU in R x R 
that are close to the traveling wave solution 

(36) M(i,x) = -max(l-e-'=(^-'=*\o) 
moving with velocity c > 0. 

Proposition 11.1. (Formal linear stability of one-dimensional traveling 
waves) 

The traveling wave u given by hSb]) is formally linearly stable with respect to equation 

Remark 11.2. In higher dimensions this result is no longer true. It is well known 
that a fingering instability occurs. However the pulsation phenomenon with which 
we are concerned in the present paper appears already in dimension 1 . 

Formal proof of Proposition 

Let us consider solutions u of (|35|) satisfying 

u{t, x) = for X < s{t) , 



(37) 



u{t,x) < for X > s{t) , 

w(t, x) — > —1 as X — s{t) +O0 , 



s'{t) = -dxu{t, s{t) -f 0+) > . 

Let us remark that a simple analysis shows that we do not have a comparison 
principle for solutions of ((37)) . 

In order to analyze the stability we transform (|37p by 

v{t,y) u{t,y + s{t)) ; 
V satisfies v{t, y) = for y < and 
' v{t, 0) = , 



(38) 



v(t, y) <0 for y > , 
^'(^: y) — ^ -1 as y ^ +00 , 
dtv = dyyV + s'{t)dyV ou R X (0, -l-oo) 

[ S'it)^-dyvit,0+) . 
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We now consider for i > a perturbation of the traveling wave 

v{y) = -max(l - £-='',0) 
with velocity c > 0. In the formal expansion 

s{t) = ct + ej{t) + 0(£2) , 

V ~ V + ew + O(e^) , 
the first order terms w{t,y),"f{t) formally satisfy 
' ^(^,0) = , 

w{t, y) — > as 2/ ^ +oo , 

dtw = dyyW + cdyW + ^'(t)dyV vn (0, +oo) X (0, +oo) , 

[ i{t)^-dyW{t,Q+) . 

Let us look for solutions of the form 



w{t,y) = e^*W{y) 



i{t) 



where Re (A) > 0. We obtain 

W" + cW ~\W ^ ce-^y , 

i.e. 

W{y) = -~!le-''y + ^A±e''^y , 

where 

Mi = -c/2 ± /4 + A and Re (^y^c^/A + Xj > c/2 . 
The function W can only be bounded if = and, by W{Q) = 0, 

A 

Finally the relation 7'(t) ~ —dyw{t^O) implies 

The unique solution of this equation is A = 0. Thus we formally proved stability of 
traveling waves. 

Acknowledgment: We thank Steffen Heinze, Danielle Hilhorst, Stephan Luck- 
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